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Abstract 



Exact expectation values of the fields e"*^ in the BuUough-Dodd model are derived by- 
adopting the "reflection relations" which involve the reflection S-matrix of the Liouville 
theory, as well as special analyticity assumption. Using this result we propose explicit 
expressions for expectation values of all primary operators in the c < 1 minimal CFT per- 
turbed by the operator $12 or ^»2,i- Some results concerning the ^»i^5 perturbed minimal 
models are also presented. 
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1. Introduction 



Computation of vacuum expectation values (VEV) of local fields (or one-point cor- 
relation functions) is important problem of quantum field theory (QFT) 0. When 
applied to statistical mechanics the VEV determine "generalized susceptibilities" , i.e. lin- 
ear response of the system to external fields. More importantly, in QFT defined as a 
perturbed conformal field theory the VEV provide all information about its correlation 
functions which is not accessible through straightforward calculations in conformal per- 
turbation theory 0. Recently some progress was made in calculation of the VEV in l-fl 
dimensional integrable QFT. In explicit expression for the VEV of exponential fields in 
the sine-Gordon and sinh-Gordon models was proposed. It was found in that this ex- 
pression can be obtained as minimal solution to certain "reflection relations" which involve 
the Liouville "reflection S-matrix" , provided one assumes simple analytic properties of 
the VEV. This result for the sine-Gordon model allows one to obtain, through the quantum 
group restriction, expectation values of primary fields in c < 1 minimal CFT perturbed 
by the operator $1,3, with good agreement with numerical data 0. In this paper we 
use the "reflection relations" to obtain the VEV of exponential fields e""^ in the so called 



BuUough-Dodd model p, 0. As is known [|T0[, for special pure imaginary values of its 
coupling constant the Bullough-Dodd model admits quantum group restriction leading to 
a c < 1 minimal conformal field theories (CFT) perturbed by the operator $1,2- We use 
this relation to obtain the VEV of primary fields in these perturbed minimal CFT. 

In Sect. 2 we present some details of the derivation of the VEV in the sinh-Gordon and 
sine-Gordon models using the "refiection relations", and show how the VEV of primary 
fields in minimal CFT perturbed by $1,3 can be obtained. In Sect. 3 we extend this 
approach and find explicit expression for the VEV of the exponential fields e°''^ in the 
Bullough-Dodd model. We show that in the semi-classical limit our expression agrees with 
known results from the classical Bullough-Dodd theory. We also run some perturbative 
checks. In Sect. 4 we study the minimal CFT perturbed by the operator $1,2- Using our 
result for the Bullough-Dodd model we propose exact formula for the VEV of all primary 
fields in these perturbed theories. We also compare our results with numerical data 
available in literature. In Sects. 5 and 6 some results and conjectures concerning minimal 
models perturbed by the operators $1,5 and $2,1 are presented. 



2. Reflection relations in the sinh-Gordon model 



The sinh-Gordon model is defined by the Euclidean action 

AshG = j d'x[ -^{d^^f + + /"e-^^ } . (2.1) 
We are interested in the expectation values of the exponential fields, 

G'(a) = (e"^),;,G- (2.2) 
As was observed in these expectation values satisfy the "reflection relation" 

G{a) = R{a) G{Q - a) , (2.3) 

G{-a) = R{a)G{-Q + a), (2.4) 

where 

Q = h-^+h (2.5) 
and the function R is related to the so called Liouville reflection amplitude S [0 , 

^2 > ^ ' \T{i-b^)J r{i-2iP/b)r{i-2iPb) ^ ' 



Note that the second of the relations (|2.4| ) follows from the flrst one if one takes into 
account an obvious symmetry of (|2.2| ), 

G{a) = G{-a) . 

No rigorous proof of the reflection relations ( |2.3| ), ( |2.4| ) is known to us. Here we give simple 
intuitive argument in support of these relations. 

Let us note that the sinh-Gordon theory (|2.1|) can be interpreted as the perturbed 
Liouville QFT in two different ways. First, one could take the flrst two terms in the 
action (|2.1|) as the action Al of the Liouville theory (in a flat 2D background metric) and 
treat the last term containing e~^'^ as the perturbation. Then naively one could write 
down the conformal perturbation theory series (expansion in the perturbation term) for 
the one-point function of ( |2.1|) , 



e''^{x))sHG = Z-'f^^-^ [ rfV.-.rfV(e"^(a;)e-^^(yi)...e-''^(2/.))L , (2.7) 

n=U 



where {■ ■ ■)l are the expectation values over the LioviUe theory Al, and Z is its partition 
function which does not depend on a. With this expression the first "refiection relation" 
( p.3|) follows from the refiection property of the Liouville correlation functions (see @] for 
the details), 

( e'^^ix) ...)l = Ria) ( e(^-")^(a;) (2.8) 

where dots stand for any local insertions. The coefficient function R{a) is related to the 
Liouville two-point correlation function 

(e"^(a:)e"'^(a:'))L= [SQ-a,a' + R{a) 5a,a'] |x - xT'"^"^""^ ; (2.9) 

its explicit form is given by ( |2.6| ). The function S in ( p.6|) can be interpreted as the 
amplitude of scattering off the "Liouville wall", as explained in [0. Alternatively, one 
could interpret the second term in ( |2.1D as the perturbation of the Liouville CFT defined 
by the first and the third terms in ( p.l| ) . Then writing down corresponding naive conformal 
perturbation theory series analogous to ( |2.7|) one would arrive at the second relation ( |2.4|) . 
In both cases the problem is that the integrals in ( p.7|) (as well as the integrals appearing 
with the second interpretation) are highly infrared divergent and therefore the naive series 
( ^TtI) does not give a viable definition of the one-point function. 

One can get around the above infrared problem as follows. Consider 2D "world sheet" 
Eg, topologically a sphere, equipped with the metric g^aix), and define a version of the 
sinh- Gordon theory on Eg with the following non-minimal coupling to the background 
metric g, 

■^Lg = I d'x^ { ^ g^-d.^d^^ + 1^ ^ + ^ : e^^ :g +^ : e"''^ :g } , (2.10) 

where Q is given by ( |2.5| ), R denotes the scalar curvature of g and the symbol : e^'"^ -.g 
signifies that these exponential fields are renormalized with respect to the background 
metric g. The first three terms in (|2.10| ) define conformaly invariant Liouville theory 
on Eg so that ( p.lO|) agrees with the first of the above interpretations of the sinh-Gordon 
model as the perturbed Liouville theory Al- Precisely this was our reason for adding the 
curvature term in ( p.lO|) . Due to its conformal invariance the Liouville theory is insensitive 
to a choice of the background metric g. If one picks a conformal coordinates on Eg, so 
that 

gua{x) = p{x) , (2.11) 
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the dependence on p{x) can be expelled from the Liouville part of the action ( p.lO|) by the 
shift 

ip{x)^ipix)-Qlogpix) . (2.12) 



This transformation brings (|2.10| ) to the form (up to field independent constant) 



^IhG = Jd'^{^ (d.^f + ^e'^ + pp'+^'^e-'^ }+Q^oo, (2.13) 

where now the exponential fields e^'"^ =: e^'"^ -.gio) are normalized with respect to the fiat 
metric gj^J = 5,^^, so that 



■g 



The term with (poo = ^iT^\x\^oo ^{x) plays no role in the perturbed theory. To be definite, 
let us take the metric g to be that of a sphere with area A, 

p(a;) = (l + 7r|x|V^)"^ ■ (2.14) 

For finite A the conformal perturbation theory for (e"''')^^^, in (|2.13|) (the expansion in 
e~'"^) makes much better sense because now the integrals analogous to those in ( |2.7| ) 
contain the factors 11^=1 [piVk) ]^^^^ providing an efficient infrared cutoff. As the result 
these calculations produce a power series of the form 



>"")f.G = ^-^^"^"-^^ E [pA'^'^r Gnia) . (2.15) 

n=2 



Owing to the property ( |2.^ ) of the Liouville correlation function each term in ( 2.15] ) satisfies 



the refiection relation ( J2.3D . Assuming that the series ( |2.15D defines a function G{a,t) 

oo 
n=2 



Yl'^=2 Gn{a) with the asymptotic 



G{a,t) ^ G{a) t^^^~-^^ as t ^ oo (2.16) 



and taking the limit A ^ oo (which brings (|2.10| ) back to (|2.1j )) one arrives at (^ 
Similarly, starting with the action which differs from ( |2.10p only in the sign of the curvature 
term, one can repeat the above arguments, this time taking the term with e~^'^ as the 
perturbation. This leads to (|2.4| ). 



Of course, these arguments do not give a rigorous proof of the refiection relation 
I) because the convergence of the conformal perturbation theory ( |2.15|) is not at all 



obvious even at finite A. Also, the existence of appropriate limiting behavior t — oo in 
( pj.l6| ) is at least problematic. On the other hand if one assumes these arguments valid the 
reflection relations ( |2.3| ), ( |2.4|) can be taken as the starting point in deriving the expectation 
values (|2.2| ). In fact, if nothing is said about analytic properties of the function G{a) the 
equations ( |2.3|) , (|2.4| ) are not nearly sufficient to determine it. However, if one makes 
additional assumption that G{a) is a meromorphic function of a, the following "minimal 
solution" to the equations ( p.3|) , ( p.4|) is readily derived 

sinh^(2a6t) 



exp 



°° dt 

T 



2sinh(62t) sinh(t) cosh ((1 + b'^)t) 



where [11 



m 



4^^ 



2 + 25^ 



;2.17) 



(2.18) 



.2+2P 2+2^2 , 

is the particle mass of the sinh-Gordon model. Note that (|2.17| ) is exactly the expres- 
sion conjectured in At the moment we have absolutely no clue on how to justify 
this analyticity assumption. We can only make a remark that while the above arguments 
leading to the reflection relations (^73|), (|2.4[) do not seem to depend on the integrability 
of the sinh-Gordon model, the simple analytic properties assumed above most likely doEI. 
However, we consider various perturbative checks of ( p.l7| ) performed in as a strong 
evidence supporting both the above arguments about the reflection relations and the ana- 
lyticity assumption. Additional support is provided by the results in Ref. where these 
assumptions are used to derive the expectation values of the boundary operators in bound- 
ary sine-Gordon model with zero bulk mass. Furthermore, in Sect. 3 we will use the same 
assumptions to obtain the expectation values of exponential fields in the BuUough-Dodd 
model. 

As mentioned in 0] , the expression (|2.17| ) can be used to obtain the expectation values 
( ^i^k ) of primary fields with conformal dimensions 

{p'l — pkY — [p' — pY 



4pp' 



(2.19) 



^ This can be compared with the situation in lattice models of statistical mechanics. While 
so called "inversion relations" (see |jl^) for the partition function can be written down for many 
models including non-integrable ones, only integrable lattice models provide enough analyticity 
for making the inversion relations a powerful tool of computing the partition functions. 
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in perturbed "minimal models" [13 

A^p/p' + A j rf\$i,3(x) . (2.20) 

This is possible because the perturbed minimal models ( |2.2C1| ) can be understood in terms 
of "quantum group restriction" of the sine-Gordon model [0, |T^, |T^, |]T^ 

Ag- J d''x\^^{d,ipf -2ficos/3ip] . (2.21) 

As is known the sine-Gordon theory in infinite space-time exhibits a symmetry with respect 
to affine quantum group Uq{sl2) with the level equal to zero and 

q = el^ . (2.22) 

The soliton-antisoliton doublet transforms as two-dimensional irreducible representa- 
tion while the bound states are scalars. The S-matrix commutes with the generators 
E± , H± , F± which satisfy the relations 

[E+,F+] = ^- [E.,F.] = ^- H+ + H. = . (2.23) 

q — q ^ q — q 



The operator Hj^ is identified with the soliton charge. Important observation made in |TG 
is that special exponential fields 

ViM^) = e*^^^(x) ( A; = 1, 2, ... ) (2.24) 

commute with the generators of the subalgebra Uq{sl2)+ G Uq{sl2), 

yi,fc(x)] = yi,fc(x)] = yi,fc(x)] = . (2.25) 

This subalgebra plays the central role in the relation between ( |2.21| ) and ( p.20|) . The space 
of states TisG of the sine-Gordon model admits special inner product [jl^ (different from 
the standard sine-Gordon scalar product) such that eI^ = (the standard scalar product 
implies = -F-). If g is a root of 1, i.e. 

= 4 , (2.26) 
P 

where p,p' are relatively prime integers such that p' > p > 1, one can isolate the 
subspace Hp G HsG consisting of the representations of Uq{sl2)+ with the spins j = 



0, 1/2, 1, . . . , p/2 — 1. The space InvlTip) of invariant tensors of Tip is identified with 
the space of states of the perturbed minimal model (|2.20|) . The relation between the 
sine-Gordon parameter // and the coupling constant A in (|2.20| ) is found in [1TT| , 



A 



(1-2/32) (3/32-1) 



r3(l-/32)r(3/3') 

r3(/32)r(i-3/32) 



(2.27) 



The theory (|2.20| ) has p — 1 degenerate ground states |0s), s = l,2,...,p — 1 |T^ which 
can be associated with the nodes of the Dynkin diagram Ap-i. The excitations are the 
kinks interpolating between these vacua, and possibly some neutral particles interpreted 
as bound states of the kinks. According to ( |2.25| ) the operators ( p.24| ) of the sine-Gordon 
model are related in a simple way to the primary fields of ( p.20| ) 



Vi,kix)=Ni,k $i,fc(x) 



(2.28) 



where Ni^k are numerical factors which depend on the normalization of The canonical 
normalization 



^i^k{x)^i,kix')) ^ \x — x'\ ^^1''' as |x — x'l ^ 



corresponds to the choice 



(2.29) 



(2.30) 



where 



7^(a) = - 



r(/3-2 + 2a/3-i) r(/32 - 2a/3) 



r(l + /32) ; r(2 - p-^ - 2a/3-i) r(2 - /32 + 2ap) ' 

is obtained from ( |2^ ) by the substitution 

b ^ i(3 , a ^ ia , fx ^ —fx . 



(2.31) 



Notice that the relation ( |2.27| ) can be written as A = —Ni^^fj.. With the normalization 
( ^) , one finds 



(0, 1 I 0,) = (-i)^('=-i) iVi,fc g{ 



1 - k 



(2.32) 



where Q{(y) is related to G{a) in ( |2.17| ) by the same substitution ( ^.31| ). The sign factor 
in ( |2.32| ) takes into account the fact that the exponential fields ( p.24| ) with even k change 
sign when if is translated by the period of the potential term in ( p.21| ) . 



For the primary fields ^i^k with / > 1 of the restricted theory ( p.20|) the situation is 
more difficult. The exponential fields 



exp 



l-l k-1 



2(3 



/3 



for I > 1 are not invariant with respect to the algebra Uq{sl2) + - Together with certain 
nonlocal fields they form finite-dimensional representations of this algebra. The calcula- 
tions become much more involved and we did not complete them yet. However, we have a 
conjecture 



(0. I ^i,k I 0,) 



sm{f\p'l-pk\) 



M 



^r(i + f) 
2r(|) 



2A, 



Qi,3{{^ + l)l-^k) , (2.33) 



where [11 



M 



2r(|) 



ttA (1-0(2^-1) 



v/^r(i + |) L (1 + ^)^ 

is the kink mass and 



r(T^)r(iny) 



P 



P — P 



(2.34) 



(2.35) 



The function Qi^sir]) for 1 3?e r] | < ^ in ( |2.33| ) is given by the integral 

• dt [ cosh(2t) sinh(t(r7 - 1)) sinh(t(77 + 1)) (r^^ 



Qi,M = exp 



-4t 



t V 2cosh(t)sinh(tOsinh(t(l+0) 2^(^ + 1) 



and it is defined by analytic continuation outside this domain 0. In writing ( |2.33| ) we 
assumed the same canonical normalization convention for the fields as in ( |2.29|) , i.e. 



/|-4A, 



(0, I $,,fc(x)$;,fc(a:') I 0,) 
Notice that (|2.33|) automatically satisfy the relation 



as 



. 



(2.36) 



(0, I I 0,) = (0, I $ 



p—l,p' — k 



Os) 



^ The expression for ( ^i^k ) proposed in Q does not contain the first factor ( 2.33 ) which carries 
the dependence on s. However, the formula in Q is equivalent to ( |2.33| ) if the expectation values 
in Q are understood not as the matrix elements between the above ground states | Os), but 
rather as the matrix elements between certain superpositions of these states which arise in the 
limit L ^ oo from the asymptotically degenerate states of the finite-size system, with the spatial 
coordinate compactified on a circle of circumference L. We will explain this point elsewhere. 
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3. Vacuum expectation values in the Bullough-Dodd model 



The Bullough-Dodd model is defined by the action f^, |§] 

Abd = j d'x[ -^{d,^? + /we''^ + /"'e"^^ } . 



(3.1) 



There is some redundancy in having two parameters /U and ^' in ( |3.1| ) because if one shifts 
the field variable in ( |3.1| ), 

(p^if + ifo , (3.2) 



n'e 2*^°, so that only the combination ^ (^')^ invariant. 



they change as /U ^ /ue^'^° , ii' 
Nonetheless, we will keep both parameters. In fact in what follows the combination 



2v/3r(i) 



^(-'- + 6+362) r(g^3^2) 



;U7rr(l + 62) 



6+3!)^ 



6+36-^ 



(3.3) 



is proven to be useful. Note that m is invariant under the shift ( p.2|) . The model ( |3.1| ) is 



integrable and its factor izable S-matrix is described in [19|. It contains a single neutral 
particle. We will show below that the mass of this particle coincides with the parameter 
m defined in (|3.3| ). In this and the subsequent sections we use the notation 



GBD{a) = {e'''^)BD 



(3.4) 



for the expectation value in the Bullough-Dodd model, where the exponential field is 
assumed to be normalized in accordance with the following short distance operator product 
expansion. 



e"'^(a;)e"^(x') 



as 



. 



(3.5) 



Notice that \a\ and \a'\ should be sufficiently small numbers, in order for ( |3.5| ) to be a 
leading asymptotic. 

Exactly as in the case of the sinh-Gordon model in Sect. 2, the Bullough-Dodd model 
can be interpreted as the perturbed Liouville theory in two different ways, with either e^'^ 
or e~i'^ taken as the perturbing operator. Using these interpretations and repeating the 
arguments in Sect. 2 which led to (|2.3| ), (|2.4| ), one arrives at two reffection relations for 



GBoia) = R{a) Gbd{Q - a) , 
GBD{-a) = R\a) Gbd{-Q' + a) 



(3.6) 
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with 



Q 



+ b. 



Q' 



2 b 
6 + 2 



(3.7) 



The function R{a) is exactly the same as in ( |2.6| ), while R'{a) is obtained from that by 
the substitution fx ^' , b ^ 6/2. As in Sect. 2, let us now assume that Gbd{(i) is a 



meromorphic function of a. Then the following minimal solution to the equations 
immediately obtained 



IS 



BD 



' 2Vr(i-62)r(i + ^) 



^ r(i + 62)r(i-f) 



m - 



r(i + eq^) r(gTj^) 



n ah—lo 



2§v/3r(i) 



exp 



T 



sinh {{2 + b'^)t) ^'(t,a) 



sinh (3(2 + b'^)t) sinh(2t) sinh (bH) 



+ 2a^ e 



-2t 



(3.^ 



where 



^{t, a) = sinh {2abt) sinh ((4 + 6^ + 2ab)t) - sinh ((2 + 26^ - 2ab)t) + 

sinh ((2 + 6^ + 2a6)t) - sinh ((2 + 6^ - 2a6)t) - sinh ((2 - 6^ + 2a6)t) 



and m is given by 



The integral in 



is convergent if 



16^ lb 
6 4 26 2 



it should be understood in terms of analytic continuation otherwise. We propose (|3.8|) as 



exact expectation values for the BuUough-Dodd model (|3.1| ). 

Expanding ( e""^ )bd — 1 + a {ip) bd + Oi^a?), one finds the expectation value of 



86 



BD 



36 



log 



fl 



' r(i-62)r(i + ^) 



^ r(i + 62)r(i-^) 



m . 



^{^+ e+sbO r(g , 3^2) 



2^v/3r(i) 



+°° sinh ((2 + 9)t) sinh ((1 - \)t) sinh ((2 + sinh ((1 + 9)t) 



(3.9) 







sinh (3(2 + 62)t) sinh(2t) sinh (6^^) 
For a = b and for a = —6/2 the integral in ( |3.8| ) can be evaluated explicitly. 



BD 



m 



BD 



24 (2 + 62) sin(^) sin (^) 



(3.10) 



These expectation values can be used to derive the bulk specific free energy of the Bullough- 
Dodd model 

Ibd = - ,Jim — log Zbd , (3.11) 
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where V is the volume of the 2D space and Zbd is the singular part of the partition 
function associated with (pJ]). Obviously, 



BD 



BD 



(3.12) 



This leads to the following result. 

fBD = 



m 



16V3 sin (^) sin (^) " 



(3.13) 



On the other hand exact expression for the specific free energy in terms of the physical 
particle mass can be obtained from exact form- factors [l2T|| , or from the Thermody- 
namic Bethe Ansatz calculations following This way one obtains exactly ( |3.13| ) 
with m understood as the particle mass. This shows that ( |3.3| ) indeed gives the particle 
mass in the BuUough-Dodd model. 



Since the above derivation of (|3.8|) is based on the assumptions (in particular, we made 
strong analyticity assumption) , some checks of this result are desirable. Simple consistency 
check is based on the known fact that for 6^ = 2 the Bullough-Dodd model is equivalent to 
the sinh-Gordon model (5TT) with = 1/2. It is possible to check that (^) calculated 
at 6^ = 2 coincides with (|2.17|) calculated at 6^ = 1/2. 

Important check can be performed in the classical limit 6^ — > 0. Consider the expec- 
tation value {e^'^)BD- In the limit ^ with a fixed, (|3.8| ) gives 



log(e>)BD 



^ (^-a' logm+Uog[^)+ £ du C{uj) ] + 0(1) 



where 



C{uj) = (21og2 + 31og3) w + log 



r(^)r(i^) 



(3.14) 



(3.15) 



On the other hand, for b^ ^ the expectation value ( 6^*^(0) )bd can be calculated 
directly in terms of the action ( |3.1| ) evaluated on appropriate classical solution ^pciix) of 
the equations of motion associated with (|3.1| ). The suitable solution depends only on the 
radial coordinate r = It can be written as ^ciix) = f {4>{f) + I ^08(2^) )' where (^(r) 
is the solution to the classical Bullough-Dodd equation 



did) + r" 



dr4>= — 



(3.16) 
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which satisfies the foUowing asymptotic conditions 

(f){r) —2a log(mr) + C(a) as r — > , 

(p[r) —>■ sin [—) sm ( — (1 — cjj Ko[mr) as r ^ +oo , 

TT o o 



where Ko{t) is the MacDonald function. The constant term C{a) in ( p.l7| ) is not arbitrary 



but must be consistently determined from the equation ( |3.16| ). Exact result for C{a) found 
in (see also [^) shows that it coincides with C(cr) defined by ( p.l5|) . Then calculation 
of the classical action gives the result identical to ( p.l4|) . 

One can go beyond the classical limit and consider the loop expansion for the expecta- 



tion values (|3.4|) . The simplest thing to study is the expectation value {(p)bd- According 



to (|3.9| ) , this quantity can be written as a power series in 6, 



( )bd = ^ log ) + 6 ( 7 + log(m/2) ) + ( 3v^7r + lOn' - 15 V' (2/3) ) + 0(6^) , 

(3.18) 

where ilj'{t) = dflog (r(t)) and 7 = 0.577216... is Euler's constant. On the other hand, it 
is possible to calculate this expectation value within the standard Feynman perturbation 
theory for the action (|3.1| ). The first term in ( |3.18|) coincides with the classical value of <f>, 
and the diagrams contributing to the next two orders are shown in Fig.l. Calculations are 



straightforward and the result is in exact agreement with (|3.18| ). 



4. Expectation values of primary fields in the minimal models perturbed by 
the operator $12 



The expectation value ( p78|) proposed in the previous section allows one to obtain the 
the expectation values 

6^cBD(«) = (e^"^)cBD (4.1) 

in the so called "complex BuUough-Dodd model" , 

AcBD = I d'x{ ^{d^? - ^e^^^ - f^'e-'^"" } , (4.2) 
which is obtained from (^T|) by replacing the parameter b by pure imaginary value b — > 



i/3 and /j — /x, n' —n'. The action (^?2|) is complex and it is not clear exactly 



how it defines a quantum field theory. Nonetheless, some formal manipulations can be 
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>f- 



>f- 



D 




-O 



>f- 




o 





Fig.l. Feynman diagrams contributing to {^)bd to the orders b and 6'^. 



done. In particular, the model ( [4 .21 ) is shown to be integrable and its factorizable S- 
matrix is identified in |T^. Although the model ( ^.21 ) is very different from ( |3.1D in its 
physical content (the model ( |4.2| ) contains solitons), there are good reasons to believe 
that the expectation values ( [4.1j ) are obtained from the expectation values by simple 
substitution 

QcBd{oi) = GBDii<y)\f,^ip ■ (4.3) 

The arguments are much the same as those which lead us to the relation between the 
expectation values in ( |2.1| ) and ( |2.21| ). Namely, it is easy to check that for fixed a the 
expression ( |3.8| ) can be expanded into a power series in b with a finite radius of conver- 
gence. This suggests that in principle ( |3.8D can be calculated by summing up the Feynman 
perturbation theory series for ( |3.1| ). At the same time the perturbation theory for ( |3.1j ) 
agrees with that for ( ^.2| ) to all orders if one makes the substitution b — > i(3. 
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The complex BuUough-Dodd model admits the quantum group restriction similar to 
the one in the sine-Gordon model []T^. There are some differences, though. The model 
( [4.2|) (in the infinite space) has a symmetry with respect to the affine quantum group 
algebra Uq{A2 ) where q is given by the same expression ( p.22| ). This algebra contains 
a subalgebra Uq{sl2) which can be used for the quantum group restriction of (0)1. If 
0^ takes the rational values ( |2.26|) the restricted theory coincides with perturbed minimal 
model [p!0| 

A^p/p' + A j (fx^i,2{x) ■ (4.4) 

with 



(2/32-1)2 r2(/32)r(l-2/32) ■ 

As in (|2.21| ), the generators E+, of the subalgebra Uq{sl2) commute with the ex- 

ponential fields ( |2.24| ) which become the primary fields ^i^k [k = l,2,...,p' — 1) in the 
restricted theory (|4.4[) . Therefore, the same arguments as in Sect. 2 lead to the following 
expression for the expectation values of these primary fields in (^ 



(0, I I 0,) = (-l)«('=-i) Ni,k GcBDi^-^P) , (4.5) 



which is similar to ( |2.32| ), except that now GcBd{o() stands for the expectation value 
in the complex BuUough-Dodd model. Here again, the canonical normalization ( |2.36| ) of 
the primary fields ^i^k is assumed. In ( |4.5|) s is an integer which labels the ground states of 
the perturbed theory. Precisely which values it takes relates to the question of the vacuum 
structure of the perturbed theory ([4. 41); we will discuss this point a little later. For now. 



it suffices to note that only the signs of the the expectation values ( [1.5|) depend on the 
choice of the vacuum. 

Particular case of (^75|) is the expectation value of the perturbing operator. 



g-2 



2WK + i)r^(i) r(|a)r(^) ^ 



,2 



2g+2 
3?+6 



(4.6) 



The above affine quantum algebra contains another subalgebra, t/ 4(5^2), which can be used 



for another quantum group restriction of (4^). We briefly discuss this case in Sect. 5. 
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where ^ is given by ( p.35| ). As usual, this expectation value is related to the specific free 
energy /i 2 of the perturbed theory ( |4.4D 



dxh,2 = {Os I $1,2 I 0,) . 



(4.7) 



On the other hand, the specific free energy for ( |4.4| ) is known exactly in terms of the mass 
M of the lightest kink present in this theory , 



/ 



sin {^^) 



1,2 



4 73 sin ( 



7r(2g+2) 
3?+6 



) 



(4. 



Combining (^.fip, (|4.7| ) and ( |4.8| ), one finds the relation between the coupling constant A 
in ( [4 .41) and the mass M, 



2^V3r(i)r(3^) 



M 



S+1 

35+6 



(4.9) 



in exact agreement with |j2^]. According to ( ^?9|) , the perturbed QFT ( ^^ ) develops a 
massive spectrum for 

0<^<1 3?eA = 0; 

(4.10) 

^ > 1 , ^mX = 0. 

In what follows we will discuss the second case only. 

As in the case of $1,3 perturbation the situation with other primary fields ^i^k with 
/ > 1 in (|4.4| ) is more difficult. However, there is a natural modification of the conjecture 
( F32| ) suitable for (|1|), 



sin ( ^\p'l — pk\) 



2«v/3r(i)r(3^) 



Ql,2((? + l)^-^fc) , 

(4.11) 



where Qi,2(??) for I 3?e 77 1 < ^ (^ > 1) is given by the integral 



QiAv) = exp 



dt ^ sinh((^ + 2)t) sinh(t(r7 - 1)) sinh(t(r7 + 1)) 



X 



t V sinh(3(^ + 2)t) sinh(2(^ + l)t) sinh(^t) 
cosh (3(^ + 2)t) + cosh ((^ + 4)t) - cosh ((3^ + 4)t) + cosh (^t) + 1 



-4t 



2^(^ + 1) 
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and is defined through analytic continuation outside this domain. The integer s in 



labels the vacuum states of ([4 .41 ). To discuss the vacuum structure of (|4.4| ) let us recall 
the basic idea of quantum group restriction. The states of ([4.2| ) can be classified according 
to the representations of the above quantum algebra Uq{sl2)- If takes the rational 



value (|2.26| ) the subspace Tip G HcBD consisting of the representations with the spins 



j = 0, l/2,l,...,p/2 — l,is closed with respect to the dynamics of ( [4.2| ). However unlike the 
sine-Gordon case, the solitons of ( |4.2D transform as the three-dimensional representations of 
Uq{sl2) with the spin j = 1. Therefore in fact there are two dynamically closed subspaces, 
Hp and containing respectively half-integer or integer spins out of the above set 
j = 0, 1/2, ...,p/2 — 1. Each of the spaces Inv(Hp) and Inv(T-C~) of invariant tensors 
associated with 7i+ and H~ can be interpreted as the space of states of certain quantum 
field theory. Therefore, this quantum group restriction of ( |3.1| ) gives rise to two different 
quantum field theories. Notice that if p is odd, the spaces Inv(Hp) and Inv(H~) are 
isomorphic because of the known property of the tensor category of representations of 
Uq{sl2) with qP = ±1, and hence the corresponding field theories are equivalent. If p is 
even, these are two really different field theories. Let us recall in this connection that the 
minimal CFT M.p/pi always has Z2 symmetry which acts on the primary fields as 

-> , . (4.12) 

If p is odd, the perturbing operator $12 is odd under the transformation ( |4.12[ ), and there- 
fore changing the sign of the perturbation in ( |4.4| ) leads to an equivalent field theory - all 
its correlation functions are obtained from those of the original theory by the substitution 
( [4.12|) . On the contrast, if p is even, the operator $1^2 is invariant under (|4.12| ); in this case 



the theories ( [4.4|) with different signs of the perturbation are expected to be essentially 
different^. Therefore it is natural to identify /nf (7i+) with the space of states of (4.4) 



with A > and Inv(T-C~) with the space of states of ( ^4.4|) with A < 0. The integers s 
labeling the vacua in (|4.11|) are related to the spins j admitted into Hp as s = 2j + 1. 
We conclude that the theory ( [4.4|) with odd p has degenerate ground states indepen- 
dently on the sign of A; however, if A > these ground states are identified with the above 
vacua I Os) with even s = 2,4, — 1, while if A < these are the vacua | 0^) with odd 
s = 1, 3, — 2. If p is even and A > there are p/2 — 1 ground states identified with 



^ The exception is the case {p,p') = (4,5) where the theories with different signs of A in (^ 
are related by duality transformation [E6[, [p7|. 
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I Og) with even s = 2, 4, — 2. Finally, if p is even and A < the theory ( [4.4]) has 
p/2 ground states | Og) with even s = 1, 3, — 1. With this understanding, the formula 
( [4 .111 ) applies to all models ( [4 .41 ), both with positive and negative A. 

In the Truncated Conformal Space method [^| was adopted to obtain numerically 



the expectation values of primary fields in the perturbed theory ( |4.4| ) for some {p,p')- It 
is interesting to compare our results to these numerical data. 

The model ( ^4.4| ) with (p, p') = (3, 4) describes the Ising model at critical temperature 
with nonzero magnetic field. In this case there is only one vacuum |0) = IO2) (we assume 
that A > 0) , and ( ^4.11| ) gives 

(0 I $1^2 I 0) = -1.27758... AT5 , 
(0 I $1^3 I 0) = 2.00314... A^ . 

According to numerical calculations in this case 

(0 I $1,2 I 0),,^ = -1.277(2) A* , 
(0 I $1,3 I 0),,^ = 1.94(6) A* . 

The VEV (0 | $1^3 | 0) was also obtained from the fit of lattice data with the result 
2.02(10) aA 0. 

The case {p, p') = (4, 5) in ( |4.4|) describes Tricritical Ising model perturbed by the 
leading energy density operator e{x) of the conformal dimension Ai^2 = jq- The minimal 
model A^4/5 contains four more primary fields (besides the above field e = $1,2 and the 
identity operator), the sub-leading energy density operator e' = $1,3 with Ai^s = | (some- 
times referred to as the "vacancy operator"), two magnetic operators a = $2,2 (^2,2 = ■^) 
and a' = $2,1 (^2,1 = jq) and $1,4 (the latter does not have obvious physical interpre- 
tation). If A < the Ising symmetry a —a is spontaneously broken and there are 
two ground states | +) =| Oi) and | — ) =| O3). If A > there is a unique ground state 
I 0) =1 O2). With these values of s ( |4.11| ) gives 



(± 


$1,2 


l±) 


= 1.46840... 


(-A)^ 


(0| 


$1,2 


0) = -1.46840 


(± 


$1,3 


l±) 


= 3.70708... 


(-A)§ 


(0| 


$1,3 


0) = 3.70708... 


(± 


$2,2 


l±) 


= ±1.59427. 


.. (-A)^ 


(0 


$2,2 


1 0) = 


(± 


$2,1 


l±) 


= ±2.45205. 


.. (-A)^ 


(0 


$2,1 


1 0) = . 
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These numbers can be compared with the numerical results quoted in 



(± I $1,2 I ±)num = 1.466(6) (-A)^ 

(± I $1,3 I ±)n™ = 3.5(3) (-A)* 
(±1 $2,2 I ±)num = ±1.594(2) (-A)^ 
(±1 $1,2 I ±)„„^ = ±2.38(6) (-A)M 



(0 I $1,2 I 0), 
(0 I $1,3 I 0), 
(0 I $2,2 
(0 I $2,1 



|0), 



= -1.465(5) A^ 
= 3.4(2) A^ 

: 
. 



The VEV ( $1,3 ) was earlier estimated in the work [ p9| as 3.78 |A| 3 , with the quoted error 
of 5 - 10% 



5. Minimal models perturbed by the operators $i. 



The afhne symmetry algebra Uq{A2 ) of ( [1.2|) contains also the algebra [7^4(5/2) as a 
subalgebra. One can use it to obtain another quantum group restriction of ([4.2|). Let p,p' 
be two relatively prime integers such that 2p < p' . As is known, for 



P 

this restriction gives the perturbed minimal model 

-Mp/p' + A / rf^x$i,5 . 



(5.1) 



(5.2) 



The above condition 2p < p' (which excludes unitary models 7Wp/p_|_i) guarantees that the 



perturbation is relevant. The coupling parameter A in (|5.2|) is related to the parameters 
in (13) as 



32 TT^ IX (//') 



l\2 



_ (4 - 5/32) (1 - /32)(4 - 3/32) (2 - /32 



2 r5(i_/g)r(M!) 
r5(^)r(i-Mi) 



(5.3) 



According to the general scheme of the quantum group restriction one expects that the 
restricted theory has a particle of the mass m (possibly among other particles and kinks) 
given by ( |3.3D with }p replaced by —(3'^ and n — /U, n' — /x'. Excluding fx{fx')'^ from 
these relations, we can expresses the perturbation parameter A in terms of the physical 
mass scale m, 



m 



2V3r(i) 



4^2 A2 (1 _ 40^(1 - 20^r2(f^) r(^) r(i=f ) 



(1+^)4 r2(^)r(^)r(^) 



12(1-5) 



(5.4) 
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Here we use the notation 



P 



As it follows from 



p' — p 

QFT (|5.2|) presumably has a massive spectrum for 



1 t 3 



53m A = , 
3?e A = . 



(5.5) 



Outside this domain the physical content of the model ( |5.2| ) is particularly unclear. We 
restrict our following discussion to the domain (|5.5|) . Then, the specific free energy of 
can be obtained from (|3.13| ) by the substitution - 

,2 



4C 



1+C 



I.e. 



/: 



m 



1,5 



(5.6) 



16v/3sin(^) sin(I^) 

Using this relation and (|5.4| ) one derives the following expression for the expectation value 
of the perturbing operator in ( ^72|) , 



:$i,5) 



7r(l-40(l-20r(^) 



r(Tfe)r(^) 

i2v/3(i + o^r(3^) sin(^)sin(^) \ r(^)r(^) 

r(fci)r(^) 



(5.7) 



8g-4 

m i+« 



To obtain more general expectations values, let us note that the generators E- , H- , F- 
of Uq4,[sl2) commute with the exponential fields e'^^^'^{x), k = 1,2,... which under the 
restriction become primary fields $i^fc(x) in the perturbed theory (|5.2|) . Using the results 
of Sect. 3 and 4 we can obtain 

:i + ^)r(|E|)r(Jd±)n2A., 



(0. I I 0,) = (-i)('=-i)- 



m 



2lv/3r(i) 

where Qi,5(?7) = Q{r))/Q{l) and the function Q{r)) for | 9fte | < ^ (^ < 1) is given by the 
integral 

dt ( sinh((l - i)t) cosh(2t?/) 



QW) ^^Py^ I \^2sinh(3(l-^)t)sinh((l + 0^)sinh(2^t) 
(^cosh (3(1 - + cosh ((3 - i)t) - cosh ((1 - 3^)t) + cosh ((1 + i)t) + 1 

2^(^ + 1) 
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In (|5.8| ), again s is an integer labeling the ground states of (|5.2|) . The ground state structure 
of ( ^.21) is not completely understood and we do not discuss it here. Note that the value 
of s affects only the sign of the expectation values (|5.8| ) with even k. 

Some of the perturbed theories ( |5.2[ ) are studied in the literature Let 
us see how our result (^]^) matches the data from these references. 

1. The model ( |5.2D with {p,p') = (2,7) was studied in In this case $1^5 = $1,2- 
Note that for this value of p there are no kinks in ( [4.4|) . However, in this case (|5.4|) agrees 
exactly with the mass m = 2Msin (j^) (where M is given by ( |4.9| )) of one of the scalar 
particles of (|4.4|) . In fact, this theory contains two particles with the masses 



m 



mi = 



2cos(fL) 



, 1712 = 1^ ■ 



The formula 



gives 



X = -i 0.0785556... m^^ , 



which is in good agreement with the Thermodynamic Bethe Ansatz calculation in ||30|| , 

^num = —i 0.0785551... m/ . 



2. The model ( |5.2| ) with (p,p') = (2,9) was studied in The model contains four 
particles with the masses 



Ttt tt tt Ttt 

mi, 7712 = 2mi cos (—), ms = 2mi cos (— ) , 7714 = 4mi cos (—) cos (—) 



30' ^15' ' ^10' '30- 

Calculations based on the Truncated Conformal Space method in this model give 



A. 



-i 0.013065... . 



Identifying m = 7712 we obtain from (|5.4| ) 

X = -i 0.0130454... mf . 



3. The case {p,p') = (3, 14) in ( |5.2D was investigated in p2|. There are six particles 
with the masses 



mi = m2 , m3 = V2mi , m4 = = 2m cos (— ) , rriQ = 2\p2m cos ( — ) . 



TT 
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The relation 



A. 



-i 0.011833... 



given in |3^] is in good agreement with (|5.4|) which gives 



X = -i 0.01183265... 



provided we identify m = m^. 

It is interesting to notice that in aU the above examples the mass m is related to the 
mass nil of the lightest particle in (|5.2|) as 



m 



mi 



2sin( 



271-g ' 

3-3C' 



(5.9) 



We believe that this is a general relation for (|5.2|) which holds as long asl/5<^<5/9. 



6. Expectation values of primary fields in tlie minimal models perturbed by 
the operator $2 1 



As is well known, the minimal models A4p/p' admit yet another integrable perturbation 



(6.1) 



Mp/p' + A / d^x^2,i{^) ■ 



The theory (|6.1|) makes sense for 2p > p'; this condition guarantees that the operator $2,1 
is relevant. For 3p > 2p' the vacuum structure of ( |6.1| ) is expected to be very similar to 



that of (|4.4| ), with p' playing the role of p. Namely, if p' is odd, i.e. the perturbation is 
odd with respect to the symmetry (|4.12| ), the theory has degenerate ground states 
which we denote | Os) with s = 1, 3, ...,p' — 2 if A < and | Os) with s = 2,4, ...,p' — 1 if 
A > 0. If p' is even and A > there are p' /2 — 1 ground states | Og); s = 2, 4, ...,p' — 2, 
while if p' is even and A < 0, there are p' /2 vacua | 0^); s = 1, 3, — 1 (compare this 
with the situation in (|4.4|) discussed in Sect. 4). The excitations are the kinks interpolating 
between these vacua (and possibly some neutral particles). The natural modification of 
the conjecture ( [4.11| ) in this case is 



(0. I ^i,k I 0,) 



sm{f\p'l-pk\) 
sm{^{p' -p)) 



2.v/3r(i)r( 
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€+1 

3^-3 



) 



2A, 



Q2M+l)l-^k), (6.2) 



where the function 



Q2,i(^) = exp 



dt ( sinh((^ - l)t) sinh(t(?7 - 1)) sinh(t(r7 + 1)) 
T 



sinh(3(^ - sinh((^ + l)t) sinh(2^t) 
(I cosh (3(^ - l)t) + cosh {{i - 3)t) - cosh ((3^ - l)t) + cosh ((^ + l)t) + l)- (6.3) 



2e(e+i) 

is obtained from ( |4.11|) by the substitution ^ — — 1 — ^, and M is the hghtest kink mass. 
It is remarkable that the function Q2,i(^) coincides with 21,5(77) in ( |5T8| ) continued to the 
domain ^ > 2. For (/, k) = (2, 1) (U) gives 



(42,1 >=(-!)' 



3V3K-l)r(S±i)sm(|li) ^r(i + i) 

1: 

3e 



X 



v^r(i)r(|^) 



3g-3 

M 2? 



(6.4) 



This formula together with the expression 

sin (Zf±i)) 



2,1 



3g-3 

4v^sin(|^) 



(6.5) 



for the specific free energy of ( |6.1| ) leads to the following relation between M and A in (^ 



2*^v/3r(i)r(^) 



M = 



TT 



TT 



A^r^(2i|i)r(i-i) 
rn%i)r(i + i) 



3^-3 



(6.6) 



which reproduces the result of p^ . 

Some checks of ( |6.2D can be made. For (p, p') = (3, 4) the model ( |6.1[ ) is just the 
off-critical Ising field theory with zero magnetic field and $1,2 coincides with the order 
parameter a. In this case all expectation values (|6.2|) agree with known result [|33|, 



(Oi 



Oi) = -(O3 



(6.7) 



r I O3) = Ms 2^ e"8 ^2 , (O2 I a I O2) = , 

where A = 1.282427... is Glaisher's constant and M = —2nX. 

In the case {p,p') = (4, 5), (|6.1| ) describes the tricritical Ising model with sub-leading 
magnetic perturbation. The theory has two degenerate ground states | O2) and | O4) (we 
assume that A > 0) and the formula (|6.2| ) gives 

(O4 I $2,2 I O4) = -1.79745... AT5 , (O2 I $2,2 I O2) = 0.68656... A^ , 
(O4 I $1,2 I O4) = 2.04451... A^ , (O2 I $1,2 I O2) = -0.78093... A* . 



(6.^ 
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In the following numerical estimates for the expectation values of these fields were 
obtained 

($2,2)num ^ -1-09 , 



(6-9) 

{^l,2)num ^ 1-2 A45 . 



Direct comparison between ( |6.8|) and (|6.9| ) is problematic because it is not clear from 



precisely which ground state is taken in calculating the expectation values (|6.9| ). The cal- 
culations in are based on Truncated Conformal Space method where one starts with 
the finite-size system with the spatial coordinate compactified on a circle of circumference 
L; the estimates are then obtained by extrapolating the finite-size results to L = oo. Let 
us note in this connection that in the finite-size system the above two ground states | O2) 
and I O4) appear in the limit L ^ 00 as particular superpositions of two asymptotically 
degenerate states (with the energy splitting ~ e~^^ where M is the kink mass) which we 
denote | /) and | //) (the first being the true ground state at finite L). Simple calculation 



which takes into account the known kink picture of the excitations in this theory (see flC 
gives the following exact relation between these states at L = 00 



2cos(fL) cos(fL) 
^ 2cos(^) , ^ cos(^) , ^ 



(6.10) 



It is easy to check using that 

(/ I $2,2 I /) = I $1,2 I /) = . (6.11) 

This result is not at all surprising as at finite L the ground state | /) is obtained by 
perturbing the conformal ground state (i.e. the primary state | $1,1)) by the operator 
$2,1- Therefore at any finite L the state | /) is a superposition of the states out of the 
conformal families [$1^1], [$2,1] and [$31]. By the conformal fusion rules all these 
states produce vanishing expectation values of both $1^2 and $2,2- In fact, the result 
( |6.11| ) can be considered as a nontrivial consistency check of our conjecture ( |6.2| ) . On the 
other hand, the expectation values 



(// I $2,2 I II) = -{II I $2,2 I I) = 1.11089... AT5 
(// I $1,2 I //) = (// I $1,2 I /) = 1.26358... A* , 



(6.12) 
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do not vanish and actually the numerical results (|6.9D appear to be reasonably close to 
(// I $2,2 I I) and (// I $1^2 I I)- We believe therefore that it is these expectation values 
that are quoted in 0. 
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